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Abstract. 

Approximations are commonly employed to find approximate solutions to 
the Einstein equations. These solutions, however, are usually only valid in 
some specific spacetime region. A global solution can be constructed by 
gluing approximate solutions together, but this procedure is difficult because 
discontinuities can arise, leading to large violations of the Einstein equations. 
In this paper, we provide an attempt to formalize this gluing scheme by 
studying transition functions that join approximate analytic solutions together. 
In particular, we propose certain sufficient conditions on these functions and 
proof that these conditions guarantee that the joined solution still satisfies the 
Einstein equations analytically to the same order as the approximate ones. An 
example is also provided for a binary system of non-spinning black holes, where 
the approximate solutions are taken to be given by a post-Newtonian expansion 
and a perturbed Schwarzschild solution. For this specific case, we show that if 
the transition functions satisfy the proposed conditions, the joined solution does 
not contain any violations to the Einstein equations larger than those already 
inherent in the approximations. We further show that if these functions violate the 
proposed conditions, then the matter content of the spacetime is modified by the 
introduction of a matter-shell, whose stress-energy tensor depends on derivatives 
of these functions. 
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1. Introduction 

The detection and modeling of gravitational radiation is currently one of the primary 
driving forces in classical general relativity due to the advent of ground-based [HEIEIS] 
and space-born detectors 5, 6^. This radiation is generated in highly dynamical 
spacetimes, whose exact metric has not yet been found. Approximate solutions to 
the Einstein equations, both analytical (such as post-Newtonian solutions [7]) and 
numerical, have served to provide insight on the dynamics of such spacetimes and the 
character of the radiation produced. These approximations, however, have inherent 
uncontrolled remainders, or errors, due to truncation of certain higher order terms 
in the analytical case, or discretization error in the numerical case. An approximate 
global solution, then, could be constructed by joining several approximate solutions 
together in some overlap region [H [9] . 

Before proceeding, we must distinguish between two different kinds of joined 
solutions: mixed ones, where one joins an analytical solution to a numerical one; and 
pure ones, where one joins two analytical solutions that have different but overlapping 
regions of validity. A special kind of mixed joined solutions have been created 
in the context of the effective-one-body formalism [TO] [11], with the motivation of 
providing accurate waveform templates to gravitational wave interferometers. Pure 
joined solutions have been discussed in the context of asymptotic matching, where 
one of the motivations is to use the joined solution as initial data for relativistic 
simulations [El OH- In this paper, we concentrate on pure joined solutions, although 
the methods and conditions we find can be straightforwardly extended to mixed joined 
solutions. 

The construction of pure joined solutions is not always simple because 
approximation methods tend to break down in highly dynamical spacetimes. The 
main difficulty lies in that approximate solutions usually depend on the existence of 
a background about which to perturb the solution. However, in highly dynamical 
spacetimes, such a background cannot usually be constructed. Thus, in those 
scenarios, the region of validity of the different approximations tends not to overlap. 
In slightly less dynamical cases, the regions of validity can overlap, but the different 
approximate metrics usually describe the spacetime in different coordinates and 
parameters. Patching |12] could be used to relate the different approximate solutions, 
but this method usually leads to an overdetermined system if we require the patched 
metric to be differentiable at the junction. A better alternative is to relate the 
approximate solutions via asymptotic matching, which guarantees that adjacent 
metrics and their derivatives be asymptotic to each other in some overlap region. 

Asymptotic matching was developed as a technique of multiple-scale analysis to 
solve non-linear partial differential equations [121,113] . In general relativity, this method 
was first studied in [HI [T5] (TB] [17\ [18] and it has recently had important applications 
to post-Newtonian theory [7], black hole perturbation theory [19] and initial data 
constructions [H [9] . Asymptotic matching requires that we compare the asymptotic 
expansions of the approximate solutions inside the region of the manifold where the 
regions of validity of the approximations overlap (the buffer zone.) This method then 
provides a coordinate and parameter transformation that relates the approximate 
metrics, such that adjacent metrics and their derivatives become asymptotic to each 
other inside the buffer zone. In essence, asymptotically matched approximations are 
guaranteed to represent the same metric components in the same coordinate system 
inside the buffer zone. 
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After matching has been carried out, there is still freedom as to how to join the 
matched approximate solutions together. The simplest way to do so is through a 
weighted linear combination of approximate solutions with transition functions. This 
method was first considered in depth in [5J [3] , where a binary black hole system was 
studied to construct initial data for numerical simulations (the so-called Frankenstein 
approach). In that work, only broad comments were made as to the type of allowed 
transition functions, requiring only that the functions be "differentiable enough," and 
the the properties and conditions these functions must satisfy were not studied. A 
priori, it might not be clear which functions are allowed such that the global metric 
still approximately solves the Einstein equations. For instance, it might seem natural 
to use Heaviside functions to join the metrics together at a hypersurface, as is done in 
the standard junction conditions of general relativity [20| . [2l l [2 ^. [23 1 f24 | [^5 ], [26 1 f2T l [28] . 
These joining procedure works well when dealing with exact solutions to the Einstein 
equations, in the sense that the joined solution is itself also a solution. As we show 
in this paper, however, when working with approximate solutions better transition 
functions need to be found. 

The purpose of this paper is to study whether a global approximate solution to the 
Einstein equations, be it pure or mixed, can be constructed directly from a weighted 
linear superposition of approximate solutions with transition functions. We thus refine, 
prove and verify many of the broad statements made in [5J [5] regarding transition 
functions. This goal is achieved by constraining the family of allowed transition 
functions for pure joined solutions via certain sufficient asymptotic differentiability 
conditions. These conditions are independent of the perturbative order to which 
asymptotic matching is carried out and the location in which the transition occurs, as 
long as it is inside the buffer zone. We then derive and prove theorems that guarantee 
that pure joined solutions constructed with this restricted family of transition functions 
satisfy the Einstein equations to the same order as the approximations. Moreover, 
we also show that if the transition functions do not satisfy the proposed conditions, 
their derivatives modify the energy-matter content of the spacetime by introducing 
a non-negligible stress-energy tensor. Finally, we extend these theorems to solutions 
projected onto spatial hypersurfaces, so that they can be directly applied to initial data 
construction schemes. These theorems then allow for the systematic construction of 
pure joined solutions and they can be straightforwardly extended to mixed joined 
solutions. 

An example of the proposed theorems and allowed transition functions is then 
provided by studying a binary system of non-spinning black holes, where the 
approximate solutions are taken to be given by a post-Newtonian expansion and a 
linearly perturbed Schwarzschild solution. We shall not perform a systematic study 
of transition functions here, but instead we pick functions that are variations of 
those chosen in [8j [9] in order to illustrate how the gluing procedure works and 
how it breaks down. We explicitly show that if the transition functions satisfy 
the proposed conditions, the 4-Ricci scalar calculated with the pure joined solution 
vanishes to the same order as the uncontrolled remainders in the approximations. In 
particular, we explicitly show that derivatives of the pure joined approximate solution 
built with appropriate transition functions are equal to derivatives of both original 
approximate solutions up to the uncontrolled remainders in the approximations. We 
also numerically show that if the transition functions violate the proposed conditions, 
their derivatives modify the matter content of the spacetime by introducing a shell 
of matter. In this manner, we explicitly verify, both analytically and numerically, 
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that a pure joined solution does represent the same spacetime as that of the original 
approximate solutions in their respective regions of validity up to the accuracy of the 
approximations used. 

This paper is divided as follows: in section [5] we review the standard junction 
conditions at a hypersurface, so that we can extend them to the case where the 
solutions used are approximate instead of exact, provided the existence of a buffer zone; 
in section [3] we study how to join approximate solutions with transition function and 
derive conditions such that a pure joined solution satisfies the Einstein equations to 
the same order as the approximations used; in this section, we also study projections of 
these pure joined solutions to a Cauchy hypersurface in order to develop conditions for 
transition functions that can be used in initial data construction schemes; in section [4] 
we study an example of the theorems formulated by considering a binary system of non- 
spinning black hole, constructing approximate global metrics with different transition 
function, and explicitly calculating the 4-Ricci scalar; in section [S] we conclude and 
point to future research. 

The notation of this paper is as follows: Greek indices range over all spacetime 
indices, while Latin indices range only over spatial indices; the symbol O(A) stands 
for terms of order A at most, while the symbol 0(A, B) stands for remainders of order 
A or B at most, where A and B are dimensionless; a tilde superscript A stands for the 
asymptotic expansion of A as defined in [T21[T5] ; the relation A = B+0(C) means that 
A is asymptotic to B with uncontrolled remainders of order C (the so-called Landau 
or asymptotic notation); we use units where c = G = 1. Symbolic calculations are 
performed with either Mathematica or Maple. 

2. Junction Conditions at a Hypersurface 

In this section, we review a variation of the standard junction conditions of general 
relativity. These conditions have been discussed extensively in the literature (see 
(SQliaiEHESlEllESlESlEIlEgand references therein.) Here, we review only those 
concepts important to the understanding of this paper, following in particular |28j . 
Certain departures from the notation of [28] are so that the generalization to the case 
of approximate solutions in the next section becomes easier. 

Let us first set up the problem. Consider a spatial (or timelike) hypersurface J 
that divides spacetime into two regions: C\ and Each of these regions possesses an 

associated metric and coordinates, {gjiJ , x< n)} an d {9$ , 2^2)}' sucn that these metrics 
solve the Einstein equations exactly in their respective regions. In the literature, this 
exact solutions and coordinate systems are also sometimes referred to as {<?^ M , £±}- Let 
us further assume that an overlapping coordinate system x a exists in a neighbourhood 
of J . The problem is to formulate junction conditions on J that guarantee that the 
joined 4-metric satisfies the Einstein equations. 

Let us make these statements more precise by considering a congruence of 
geodesies piercing J defined with respect to first metric in C\ and the second metric 
in Ci (see [28] for a detailed definition of such congruence.) Let then I denote the 
proper time (or proper distance) along the geodesies, such that I = corresponds to 
when the geodesies reach J. The joined solution then takes the following form [28] : 

v-e(%^ + 0H)ff^ (1) 

where O(-) is the Heaviside function. Equation ([!} implicitly uses coordinates x a that 
overlap both the coordinates local to C\ and C 2 in a neighbourhood of J. 
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We now proceed to formulate the junction conditions. The first junction condition 
arises by requiring that the 4-metric be continuous across J ', in a coordinate system 
x a that overlaps both x" and x^ in an open region that contains this hypersurface. 
This condition can be expressed in a coordinate-independent way by projecting it to 
J . Then, in terms of 3-tensors this condition becomes 



h 



j 



.,(2) 
'ab 



J 



(2) 



where here h a b is the 3-metric associated with the junction hypersurface, i. e. 

hfiv = g^ u + e±n fl n u , (3) 

with n a normal to J and e± — ±1 depending on whether J is spatial (— ) or timclikc 
(+.) Equation guarantees that the hypersurface J has a well-defined geometry. 
These equations also imply that the metric is differentiable across J, except for its 
normal derivative to J that in general is discontinuous. 

The second junction condition arises by requiring that this normal derivative does 
not lead to violations of the Einstein equations across J . In terms of 3-tensors on J7 - , 
this condition becomes 



K 



J 



= K 



J 



(4) 



One can show that the failure of these equations to be satisfied changes the distribution 
of energy-momentum tensor of the spacetime and gives raise to a shell of matter with 
stress-energy tensor [25J 



S a b — 7T" 
S7T 



K 



J 



- K 



(2) 



j 



J 



(5) 



In the next section, we shall be mostly interested in vacuum spacetimes, for which 
such a stress-energy tensor should vanish. 

The satisfaction of the Einstein equations out of J by JTj and the absence of a 
stress-energy tensor as in ([5]) then guarantees that the junction conditions and 
([4|] are also satisfied. Exact solutions, however, are rarely available for astrophysically 
realistic scenarios. In that case, one must rely on approximate solutions, for which 
similar conditions to those discussed here can be found, as we shall study in the next 
section. 



3. Pure joined solutions 

In this section we build a pure joined solution by extending the standard junction 

.(12) 

conditions to the case where the metrics g^u are only approximate solutions to the 
Einstein equations. For simplicity, we assume a vacuum spacetime and that there 

(1 2) 

exists analytic approximate expressions for g^J , such that pure joined solutions are 
sought. The conclusions of this section, however, can straightforwardly be extended 
to other cases, where numerical solutions are available instead of analytical ones, 
provided information about first and second derivatives of the numerical solutions is 
also available (note that for numerical solutions the continuum derivative operator 
must be replaced by its finite counterpart.) As shown in J5J[5], the first step in joining 
approximate solutions is to apply asymptotic matching inside some overlap region. 
Once this has been done, one can search for conditions such that the pure joined 
metric tensor satisfies the Einstein equations to the same order as . We here first 
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present the basics of asymptotic matching as applicable to this problem [5] • We 
then search for asymptotic junction conditions (or buffer zone conditions), which are 
asymptotic in the sense of [THEI] anc h thus, are to be understood only approximately 
to within some uncontrolled remainder. 

3.1. Asymptotically matched metrics 

Consider then a manifold M. that can be divided into two submanifolds with boundary 
C\ and C*2 , each equipped with a approximate metrics g$ and gj?] and a coordinate 
system x a and x a respectively. These metrics are approximate in the sense that they 
solve the Einstein equations to £)(e„ n ) for some I € N in their respective submanifolds, 
e.g. in vacuum 

G^\g$]=0(& +1 ), (6) 

where n = {1,2}, {l\,t-i) are real numbers greater than zero, and e„ -C 1 is 
some dimensionless combination of parameters and coordinates relative to the n-th 
submanifold. The symbol C(e^" +1 ) refers to terms of relative order e^" +1 with respect 

(n) 

to the leading order term in the approximate solution g^v . In principle, there could 
be logarithms of e n present in the remainders, such as in high-order post-Newtonian 
expansions, but we neglect such terms here because they shall not affect the analysis 
of this paper. Notice that we use here bars to denote the different coordinate systems 
(as opposed to numbers, as in the previous section) because we must be more careful 
about the coordinates used in these approximate metric components. For concreteness, 
let the region of validity of g$(x a ) be defined by x a » xf n and that of g${x a ) 
by x a <C x" Mt . These inequalities define a spacetime region of validity, since the 
approximate metric might not be valid for all times. For example, such is the case 
for the post-Newtonian metric of two point particles in quasi-circular orbit, which is 
valid only for times t -C i C) where t c is the time of coalescence. 

Let us further assume that these submanifolds overlap in some 4- volume, defined 
by the intersection B = C\ n C2, and sometimes referred to as the overlap region 
or buffer zone. The boundary of the buffer zone, dB, cannot be determined exactly, 
because it is inherently tied to the regions of validity of the approximate solutions, 
which themselves are only defined approximately. With this in mind, let us further 
assume that the charts {x a } and {x a } are defined in the neighbourhood of any field 
point in the buffer zone and that they satisfy x a = x a + 0(e i 1 1+1 ,e e 2 2+1 ). The buffer 
zone can then be asymptotically defined by the following condition: {x^x™) 1 / 2 <C 
{x a x a ) 1 / 2 -C (^out^o 1 *) 1 ^ 2 > where indeces are raised or lowered with the local metric 
to Ci t 2- Wherever possible, we use the Landau notation, which specifically specifies the 
behavior of the remainder. The definition of the boundary of the buffer zone should be 
understood only in an asymptotic sense, as defined in [THE]. This boundary is made 
up of two disconnected pieces, dB- and dB + , defined via x a = xf n + 0(e e 1 1+1 , e^ 1 ) 
and x a = Xg Ut + 0(e 1 1+1 , e^ 1 ) respectively. The definition of the buffer zone can be 
thought of in terms of a simplistic spherically symmetric static example, where x a — > r, 
dB is a spherical shell and dB± are 2-spheres. However, in practical examples, such as 
binary black hole spacetimes, the boundary of the buffer zone is not simply a spherical 
shell, but instead it acquires some deformation in accordance with the deformation of 
the spacetime that is being modeled. 

Before proceeding with the description of asymptotic matching, let us make 
some comments on the approach adopted in this paper. In the previous paragraphs 
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and in what shall follow, we have not adopted a rigorous geometrical approach in 
the description of asymptotic matching, the definition of the buffer zone and the 
submanifolds. A rigorous geometrical approach, for example, would impose other 
constraints on B, such that it is properly embedded in C\ and Ci- Furthermore, 
such an approach would describe the conditions under which asymptotic matching 
can be performed, since in general this method is valid only locally. Instead of such 
an approach, in this paper we adopt an analytical one, where we only provide a 
minimal amount of details in order to make the presentation simpler. For example, 
we shall restrict our attention to coordinate grids that coincide when both e\ — ► 
and 62 — > and, henceforth, we shall call regions instead of submanifolds. Such 
an approach is adopted because the aim of this paper is not to provide a detailed 
geometrical account of asymptotic matching, but instead to study transition functions 
and approximate joined solutions (see [29] and references therein for a more detailed 
and rigorous geometric account of asymptotic matching.) 

The approximate metrics live in different regions and depend on different 
coordinates (x a ,x a ) and local parameters (8 a ,9 a .) Examples of these parameters 
are the mass of the system and its velocity. However, since both g$ are valid in 
B, both coordinate systems must be valid in the buffer zone (i.e., the charts of C n 
overlap in B.) Using the uniqueness theorems of asymptotic expansions [121 113] . one 
can find a coordinate and parameter transformation to relate adjacent regions inside 

B. In order to achieve this, one must first compute the asymptotic expansions of the 

~ 2 

approximate line elements, ds(„), near the boundaries of the buffer zone, dB, and then 
compare them inside B but yet away from dB, i.e. 

dsl ) -ds 2 i2) = 0(e t 1 i +1 ,ef +1 ), in B\dB, (7) 

where the \ symbol is the standard exclusion symbol of set theory and where 
0(e l i +1 , f-2 2+1 ) stands for uncontrolled remainders of order e^ 1+1 or e l 2 2+1 as defined at 
the end of Sec. 1. Using ([7]), a coordinate transformation ip : x a — ► x a and a parameter 
transformation <fi : 9 a — > 9 a can be found inside B. We have written ([7]) in terms of the 
line element, but we could have easily written it in terms of the metric components. 
In fact, depending on the calculation, it might be necessary to asymptotically match 
different components of the metric to different order. For example, for the construction 
of initial data, the goo and gij components of the metric need to be matched to lower 
order than the goi components [9] . In what follows, we will assume that all metric 
components have been asymptotically matched to the same order. When this is not 
the case, the conditions that we shall propose on transition functions can be adapted 
by noting that one must use the highest order to which metric components have been 
asymptotically matched. 

These transformations guarantee that all components of adjacent 4-metrics are 
asymptotic to each other inside the buffer zone [T21 Q2| . This fact implies that the 
derivatives of adjacent metrics are also asymptotic to each other 

_ ° (2) _ (nrji+i-m i 2 +i- m] f8 s 

dx^...dx^ 9 ^ dx^...dx^ g ^~ (l ' 2 h {) 

for all m < I. However, the relative order to which the derivatives of adjacent metrics 
are asymptotic to each other is in general not the same as the relative order to which 
the adjacent metrics themselves are asymptotic to each other. This reduction in 
matching accuracy is due to the implicit assumption that asymptotic matching is 
carried out to some finite order. In ([8]), we have assumed that this decrease in accuracy 
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occurs in single powers of 6^2, but if this is not the case the results of this paper can 
be rescaled appropriately. 

Let us provide a general example of such matching accuracy reduction, without 
specifying a particular spacetime (for a more detailed example see [HI - ) Consider 
a post-Newtonian expression that is known to C(l/c 4 ) and that is asymptotically 
matched to a perturbed black hole solution in some buffer zone to 0(l/c 4 , e 2 ), where 
c is here the speed of light and e is the black hole perturbation theory expansion 
parameter. Let us now take the time derivative of the post-Newtonian expression, 
which results in a term of C(l/c 5 ). In post-Newtonian theory, one would never 
truncate the differentiated term at 0(l/c 4 ). However, in the theory of asymptotic 
matching to a finite order, the time derivative of the post-Newtonian quantity remains 
matched only to 0(l/c 4 ,e 2 ), and thus it will disagree with the perturbed black hole 
quantity at 0(1 /c 5 ). It is in this sense that differentiation usually decreases the 
relative order to which two expression have been asymptotically matched (we should 
note that when asymptotic matching is carried out to all orders, then this reduction 
in matching accuracy disappears [29].) 

Let us provide a more explicit example that, although similar in spirit to the one 
described above, does not require general relativity. Consider two functions 

where /o,i are constants and x(ivt) and y(£lt) are periodic funtions of time with period 
1/lo and 1/VL respectively. Let us pretend that f(t) and g{t) are approximate solutions 
to the same differential equation in the limits e± <C 1 and £2^1 respectively. If we 
asymptotically match these functions in a buffer zone where both e± ^ 1 and e2 -C 1 
up to uncontrolled remainders of 0(e 2 , e 2 ,), we discover that at t = the constants are 
/o = 1 and fi = j/(0)/[9x(0)]. If we now study their time derivatives, we find that at 
t = 

df 2/(0) (dx\ dg e 2 fdy\ 

with uncontrolled remainders of 0(e 2 , e 2 ,). Note that the derivatives are not equal to 
each other to this order because the velocities uj and have not yet been determined 
via asymptotic matching. The reduction in matching accuracy is here explicit since, 
while f(t) and g(t) have been asymptotically matched up to uncontrolled remainders of 
0(e 2 , e 2 .), their time derivatives match only up to uncontrolled remainders of C(ei, e 2 ). 

Henceforth, we assume that asymptotic matching has been carried out and that 
one of the metrics has been transformed according to {-0, <f>} such that ([7]) holds. 
We refer the reader to [T^l H3] for more details on the theory of asymptotic analysis 
and [3 [141 [T5l [16l [17l [18] for a more detailed discussion on asymptotic expansions 
in general relativity. 



3.2. Asymptotic Junction Conditions 

Let us now return to the buffer zone and note that it can be foliated by a family of 
junction hypersurfaces J7i, where i labels the member of the family. For simplicity, 
we choose these hypersurfaces to be timelike, since in the next section we project the 
4-metric to a Cauchy (spatial) hypersurface £ and it is convenient then that Ji be 
orthogonal to E. Asymptotic matching has provided a coordinate transformation to 
relate the charts inside B, so, in particular, these matched coordinates are valid in 
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an open region containing every Ji as long as this is not close to dB. We look for 
asymptotic junction conditions in this subregion of B {i.e. away from dB), in terms 

of differentiability conditions of gffl . 

Furthermore, let us also consider a family of geodesic congruences 7, that pierce 
Ji. The ith member of the family is parameterized by proper distance l\ to the ith 
hypersurface J7i, such that £i = occurs when that member reaches that hypersurface. 
Such a family of geodesic congruences is defined with respect to the approximate 
metrics gffi since these metrics are equal to each other up to uncontrolled remainders 
in the asymptotic matching scheme. Even though these geodesies are not strictly 
necessary for the construction of asymptotic junction conditions, we find them useful 
to define a measure of distance to Ji, which shall later become important in the 
definition of transition functions. 

We can now begin to look for asymptotic junction conditions by considering the 
following 4-metric tensor: 

9^=F{t i )g^ + [l-F(i i )] gfl (11) 

where F{-) is a (proper) transition function that smoothly ranges from zero to unity 
inside B and which will be defined more rigorously later. Equation (fTTj) is motivated 
by the fact that any tensor E aia2 ' " ai blb2 ... bj can always be split into 

E aia2 - a ' blb2 ... b] = G aia2 --- ai b 1 b 2 ...b J F(£ i ) 

+ [l-F(t i )]H a ^-^ blb2 ... bo , (12) 

provided that 

& b^.-.bj — Lr b 1 b 2 ...b j — tl b 1 b 2 ...b j (±0) 

inside the transition region and that the proper transition function is sufficiently 
regular. The regularity requirement is to guarantee that both terms of (|12p are 
differentiable. In (|11[) this split is valid because g\}} is asymptotic to gjtf in B up 
to uncontrolled remainders once the maps {ip, <p} have been applied. As we will see 
below, however, care must be taken when constructing such proper transition functions 
F(-) to avoid ruining the differentiability properties of the joined metric. 

Let us now define proper transition functions in order to clarify how to merge the 
metrics via (jlf p . A proper transition function is a smooth real map J- : !ft — * [0, f] 
that ranges from zero to unity inside some transition window w, while it acquires the 
value of 1/2 as li — > 0, and that satisfies the following conditions: 

as £ t ^£+, 

F(4)->0 as £i^£-, (14) 

where here £± is the proper distance to dB±. Note that the point £i = is where the 
global approximate metric of pip contains equal contributions from g£} and gffi . The 
transition window will be studied later in Sec. 4, but we can think of it qualitatively 
as a parameter of a proper transition function that determines the region where 
these functions are significantly different from unity or zero. The requirement that 
a proper transition function vanishes or tends to unity at dB± is necessary to avoid 
contamination of g\}} in Ci and vice versa, since in general gjji) has large uncontrolled 
remainders and could diverge outside its region of validity. In fact, the speed at which 
a proper transition function must tend to unity or zero will depend on the speed of the 
growth of the uncontrolled remainders of the approximations outside their region of 
validity, as we shall study later. This definition does not constrain how the transition 
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function behaves inside the transition region. Also note that this family of functions 
tends to the Heaviside function of the previous section as w — > and that it does not 
need to be analytic. We shall not constrain this family further for now, but instead 
we search for conditions on this family such that (jlip satisfies the Einstein equations 

to the same order as gjlu ■ 

With these proper transition functions, we immediately see that the joined metric 
itself also satisfies the following set of asymptotic conditions: 

9»»=9$+Oi4 1+1 A* +1 l in C,\B, 

9^=9$+0(e[i+\ei* +1 ), in C 2 \B, (15) 

while in the buffer zone the metric is some weighted linear superposition of both 
approximate solutions. Thanks to asymptotic matching, the approximate metrics are 
identical inside the buffer zone up to uncontrolled remainders [see JSJ] and, thus, this 
linear superposition is valid there, in spite of the non-linearity of the Einstein field 
equations. Also note that we are free to choose any junction hypersurface J% to join 
the metrics as long as it is inside B but away from dB. However, there is usually 
a typical choice of Ji, given by the surface where the error bars of the approximate 
metrics become comparable. Such a choice is not unique, but has previously proved 
to be close to optimal in certain scenarios [3 [9] . Asymptotic matching then seems to 
be a good technique for the construction of a pure joined 4- metric as given in (fTT|) . 
as long as we find transition functions that are sufficiently well-behaved so that their 
derivatives do not introduce errors larger than those contained in the approximations. 

We now proceed to determine the asymptotic junction conditions by analogy with 
the standard junction conditions discussed in the previous section. The first junction 
condition of the previous section is automatically satisfied asymptotically in B via 
(0 with uncontrolled remainders of 0(e e 1 1+1 , e^ 1 )- However, the second junction 
condition is not necessarily satisfied because (0 does not guarantee differentiability 
across Ji to the same order as continuity. Differentiating (fTTjl we obtain 

9^, a = F(£ t ) 5 W Q + [1 - F(£i)] gW a + F >a {g^} , (16) 

where we have defined the operation {E} = E^- 1 ' — E^ 2 ' for any function E. In order 
for pT|) to be a solution to the Einstein equations, we should require that the third 
term be as small as the uncontrolled remainders of the first two terms. By JSJ), we 
know that the last piece of the third term is of 0(e{ 1+1 , e^ 1 ). On the other hand, 
the first two terms are bounded below by their smallest size, which is of C(e^,e 2 2 ). 
We, thus, arrive at the condition 

F a = 0(e°,e°), in B. (17) 

One can show that this condition is sufficient, since it excludes cases where F a = 
0(ln(ei, e 2 )). With this condition, (fl7j|) becomes 

<w = m) g$,« + [i - nm 9% + ®(4 i+1 , 4 2+1 )- m 

In vacuum, another condition must be imposed on the 4-mctric in order for 
the Einstein equations to be asymptotically satisfied inside B. This condition 
can be enforced by requiring that the 4-metric tensor be asymptotically C 2 in B. 
Differentiating (fT6| we obtain 

g^ a p = F&) gW af3 + [1 - F(ti)] g™ af) 

+ 2F da {g^ lP) }+F al3 {g^}. (19) 
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Here the parenthesis on the indices represent the standard symmetry operation 
9(a\b,\c) — 1/2 (fjab,c + fjcb.a)- In order for the metric to be asymptotically C 2 in B, 
we must require that the last two terms be much smaller than the first two. This 
requirement can be enforced by requiring that 

F, a =0(e 1 ,e 2 ), in B, 

F a/3 = 0(e?,e°), mB. (20) 

The first condition in (|20|) is a refinement of (fT7|) . while the second condition is new. 
These equations are also compatible with F^ a <C F tCt p 1 which is a consequence of the 
fact that {g^u.a} — i^ 2 )- As one can show, the decrease in matching accuracy 

of the derivatives of matched expressions has lead to different conditions for the first 
and second derivatives of the transition functions. In this way, the last two terms of 
P|) become of 0{e[ 1+1 ,e 2 2+1 ) or smaller, because by ©, {g^} = 0(e[ 1+1 7 e e 2 2+1 ) and 
{g^v,a} = 0{e^,e^). With these conditions, (fTO|) becomes 

g^ap = F(£i) gU afl + [1 - F(l t )} gV afl + 0(£ + \4> +1 ). (21) 

Transition functions that satisfy (j2"0|) are enough to guarantee a C 2 global 
metric in B, irrespective of the matching order or the junction hypersurface chosen. 
Requiring the joined solution to be C 2 in B guarantees that the Einstein tensor is also 
continuous across any Ji. Since we are dealing with a vacuum spacctime, the junction 
hypersurfaces do not represent any physical boundary, such as a shell of matter, and 
can thus be chosen arbitrarily as long as they are inside B and away from OB, so that 
the approximations are still valid. 

We have then derived conditions on the 4-metric and, thus, on the transition 
functions, that guarantee that the Einstein tensor is continuous in B. The standard 
junction conditions in the presence of matter do not necessarily require the metric to 
be C 1 across J, whereas here we must require it to be at least asymptotically C 2 in B, 
so that no artificial features are introduced at the boundary. We can now formulate 
the following theorem: 

Theorem 1 (Buffer Zone Junction Condition). Consider a spacetime manifold M. that 

(1 2) 

can be divided into two submanifolds with boundary C\ and C<z, inside which g^ are 
approximate solutions to the vacuum Einstein equations up to uncontrolled remainders 
of 0(e l 1 1+1 , e 2 2+1 )- Let the boundaries of the submanifolds be defined asymptotically by 
the approximate boundary of the region of validity of g^u and let the intersection of 
these submanifolds, B = C\ n C 2 , be called the buffer zone. Let us foliate the buffer 
zone with a family of timelike hypersurfaces Ji and consider a family of geodesic 
congruences 7.;, such that the ith member is parametrized by proper distance £i to the 
ith hypersurface Ji. Consider the family of joined ^-metric tensor parameterized by 
ii, namely 

gVbV = F(£i) 5 « + [1 - FW)] g$, (22) 

where F(-) is a proper transition function as defined by Then, equation 

is also an ti-independent approximate solution to the Einstein equations up to 
uncontrolled remainders of 0(e e 1 1+1 , e 2 2+1 ) if and only if the following conditions hold: 

(i) The metrics gfy have been asymptotically 

matched in B up to uncontrolled remainders 
ofO(et+\ei 2 +% 
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(ii) F Q =C(e 1 ,e 2 ), in B, 
(tit) F tafi = 0(^,4), m B. 

Proof. The proof follows directly from the calculation of the Ricci tensor, which should 
vanish in vacuum. Let us then use the joined metric to calculate the Christoffel 
connection. Doing so we obtain 

r^ = r; 7 (1) F(^) + r; 7 (2) [i-F(^)] (23) 

+ \d aS (F 7 {gps} + F {g 7S } - F t s {g M }) , 

where the {•} operator was defined in P^|) . We have here used the fact that the 
inverse metric can be written as g af3 = F g^ + (1 — F) g"^, neglecting terms that are 

proportional to F(l — F) because they are of 0(e\ l1 , e 2 / 2 ). Outside the buffer zone 
the terms in parenthesis clearly vanish, but inside this region they could be large. 
However, note that by condition (i), {g^} = 0(e l ^ +1 ,e l 2 2+l ), while by condition (ii) 
the derivative of the transition function satisfies F a = 0(e\,€2)- Thus, the term in 
parenthesis satisfies 

(F n {gps} + Fp {g lS } - F 5 {g 01 }) « 0(e^+\ e| 2+1 ). (24) 

Due to the precision of the Christoffel symbols, namely that ^ < a~ = r^ 2 -* in B up 

to uncontrolled remainders of 0(e e 1 1+1 , e e 2 2+1 ), combined with conditions (?) and (ii), 
we are allowed to write the connection as 

r| 7 = r; 7 (1) F(^) + r; 7 (2) [1 - f^)] + o&+\4 +1 ). (25) 

However, in order to compute the 4-Riemann tensor, we need the derivative of 
the connection. This quantity is given by 

r ^ = r; 7 ( >(^) + r; 7 (2 J [1 - F(e,)} + f„ {r^} (26) 

+ \g a& {F^a {gps} + F t<T p {g lS } - F M {gp^}) 

+ \g a& (F n {gps,*} + Fp {g^ a } - F s {gp 7 ,v}) , 

where other terms either vanish by (|24p or are negligible. Equation (|26p must be 
obtained by differentiating ([23]) . rather than (f25|) . because otherwise we would miss 
the second derivatives of the transition functions. The third and fourth terms of ([26| 
are already of C(e^ 1+1 , ^i +1 ) by ((HJ), condition (i), (ii) and (Hi). The last term is 
finally also of 0(e{ 1+1 , e l 2 2+1 ) because although {g^v,*} = 0(ei,e 2 2 ), by condition (ii) 
F >a = 0(ei, €2). We are therefore left with 

r ^ = r; 7 ( >(^) + r;^ [1 - f^)] + o($ + \4' +1 ), (27) 

which again leads to the asymptotic condition r^J-^ — T^ 7 ^ = 0(e{ 1+1 , £ 2 2+1 ) in B. 

We are now ready to compute the Ricci tensor. Once more, using Equation (|25p 
and (j2"T|) we can write 

Rc0 = R§F(£ i )+R% [l-F(£ i )} + 0(e{^\e e i +1 ), (28) 



Frankenstein's Glue: 



Transition functions for approximate solutions 13 



where the third term groups all the cross terms that are 0{4 1+1 ,£ +1 ) in B. 
Furthermore, we know that the approximate solutions g^u satisfy the Einstein 
equations to C(e„"), which then implies that 

< } =0(4" +1 )- (29) 
We thus arrive at the conclusion that 

R a p = 0{e{- + \/i +1 ), (30) 

which then proves the theorem. □ 

This theorem allows for the construction of pure joined solutions, with a restricted 
class of transition functions. Although the theorem has been formulated for vacuum 
spacetimes, it also holds for non-vacuum scenarios as discussed above. The proof 
for the non-vacuum case can be established simply by following the above proof and 
realizing that now — 8ttT^ = 0(e^ n+1 ). Also note that the ideas of this section 
can be extended to mixed joined solutions, by replacing condition (i) by some other 
condition that guarantees that the approximate solutions represent the same spacetime 
in B. Finally, note that the conditions we impose on the transition functions are not 
very stringent, thus allowing for a wide range of possible functions. 



3.3. Projection to a Cauchy Hypersurface 

Let us now specialize the pure joined solution of (jlip to joined initial data on a Cauchy 
hypersurface. Since this hypersurface is by definition spacelike, it is convenient to have 
a foliation of the buffer zone that is timelike as given in the previous section. The data 
constructed in this section consists of an induced 3-metric on the Cauchy hypersurface 
and its extrinsic curvature. 

Consider then a Cauchy hypersurface £ E A4 , on which the 4-dimensional regions 
(Ci,C2, B) become 3-dimensional surfaces. We still have a foliation of the buffer zone 
by an infinite number of timelike junction hypersurfaces Ji of S. These hypersurfaces 
are now actually submanifolds with boundary of co-dimension 1 with respect to £ but 
co-dimension 2 with respect to M.. 

In either C n we can now define the 3-metric (h^) and the extrinsic curvature 

(K-ab ) °f ^ by projecting the asymptotically matched approximate metrics g^t) to 
this hypersurface [30 . In analogy with the previous section, let us then define these 
objects via 

h ab =F(i i )h§ + [l-F(t i )] h§, 

K ab = F{£i) K$ + [1 - F{li)\ K™. (31) 

In (|3"Tj) . F(-) is a proper transition function as defined by HH). Provided this 
transition function satisfies conditions (ii) and (Hi) of Theorem 1, then (|31[) satisfies 
the constraint equations of general relativity. 

There might be some concern that the extrinsic curvature of (|3"TT) does not 
correspond to the same hypersurface S as that described by the 3-metric because 
derivatives of the transition functions have been neglected. However, if the transition 
functions satisfy the conditions of Theorem 1, then these derivatives are of the same 
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order as the uncontrolled remainders. In order to show this, we can compute the 
extrinsic curvature from the spatial metric of (|31[) directly, i.e. 

K ab = F(li) K$ + [1 - F(£i)] - \ {h ab } CftFiii), (32) 

where is the Lie derivative along the normal vector to E. Clearly, since F(-) 
is a scalar function, the Lie derivative reduces to the directional partial derivative 
of this function along the normal vector, i.e. — n a F yCc . Here n a could be 

that associated with either of the approximate solutions, since in the buffer zone 
n Ti) ~ n ?2) = C( e i +1 ) e 2 2+1 )- Since the 3-metrics are already asymptotic to each other 
up to uncontrolled remainders of C(e^ 1+1 , e l 2 2+1 ) by ©, we must only require that the 
Lie derivative be of 0{t\, e-i). This condition is consistent with the conditions of the 
asymptotic junction theorem. 

We see then that the global 3-metric and extrinsic curvature of (j3"Tj) represent 
the same data as that obtained from (jlip directly up to the uncontrolled remainders 
in the approximations. We can then formulate the following theorem, which can be 
viewed as a corollary of Theorem 1 : 

Theorem 2. Consider a spacetime manifold M with approximate metrics g$ and 

(2) 

g)iv that satisfy the vacuum Einstein equations up to uncontrolled remainders of 
C(e^ 1+1 ) and 0{e t 2 2+1 ) on submanifolds with boundary C\ and C 2 respectively. Let 
these submanifolds intersect on a 4-volume B = C\ n C 2 and foliate the 3-dimensional 
projection ofB onto a Cauchy hyper surf ace E with timelike junction hypersurfaces Ji. 
Consider the 3-metric and extrinsic curvature of S constructed via 

h ab = F(ti) h§ + [1 - F(ii)] h%, 

K ab = F(ti) K$ + [1 - F(£i)] K$, (33) 

where F(-) is a proper transition function as defined via {1$ , &i is the proper distance 
to Ji on S, and {h^ b ' 2 \ K^ 2 * 1 } are the 3-metric and extrinsic curvature o/E associated 
with g^u" 2 ^ ■ Then, the 3-metric and extrinsic curvature (33\) satisfy the constraint 
equations of General Relativity on E to the same order as gfy 2 ^ if and only if the 
following conditions are satisfied: 

(i) The metrics gfy 2 ^ have been asymptotically 
matched in B up to uncontrolled remainders 
o/0(^+V/ +1 ), 

(ii) F )Ct = C(ei, €2), in B, 

{Hi) F^ = O(e1,e 2 ), m B. (34) 

Proof. The proof of this theorem is established by projecting the 4-metric onto E. 
The 3-metric is given by 

h a /3 = g a /3 + n a np, (35) 

where n a is the covariant normal vector to E. Using (JTTJ) and the fact that any smooth 
tensor can be decomposed with transition functions via (|T2p . we rewrite the spatial 
part of the 3-metric as 

ha b = h$F(e i ) + [l-F(e i )}h<$. (36) 
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Table 1. Description of the division of the spacetime into zones. 



Zone 


' in 


rout 




Inner zone BH 1 (Ci) 







fi/b 


Inner zone BH 2 (C 2 ) 





< b 


f2/b 


Near zone (C3) 


> m A 


< A/27T 


m A /r A 



The extrinsic curvature is given by 

K a b = £nhab- (37) 

Inserting p6[) into (|37p we can rewrite the extrinsic curvature as given in (|32p . 
However, since the transition function satisfies F a = 0(£i, £2)1 the last term of that 
equation can be neglected and we obtain 

K ab = K^F(£ t ) + [1 - F(ti)] K%\ (38) 

Now, recall that the 4-metric satisfies the vacuum Einstein equations up to 
uncontrolled remainders of 0(e^ 1+1 , e 2 2+1 ) by Theorem 1. Therefore, this metric also 
satisfies the constraint equations, since these are related to the temporal components 
of the Einstein tensor and the normal vector to E. Recall here that the normal can be 
that associated with either approximate solution, since these vectors are asymptotic 
to each other inside the buffer zone. Since the 3-metric and extrinsic curvature come 
directly from a projection of this 4-metric, it follows that this data must also satisfy 
the constraints to the same order. □ 



4. A simple example 

In this section, we investigate the impact of different transition functions on the 
satisfaction of the Einstein equations. For this purpose, we pick a metric that has 
already been matched in 9], henceforth paper I. This metric represents a binary 
system of Schwarzschild black holes in a quasicircular orbit. We here show explicitly 
that as long as the transition functions satisfy the conditions of Theorems 1 and 2, the 
joined metric satisfies the vacuum Einstein equations and, thus, also the constraint 
equations. 

As explained in paper I, the manifold can be divided into 3 regions (table [J). The 
symbols of table Q] represent the following quantities: C n labels the nth regions, where 
n = {1, 2, 3}; A labels the black hole, with A = {1, 2}; ri n and r out are the approximate 
inner and outer boundary radius of each C n regions as measured from the Ath black 
hole; and are the radial distances as measured from the Aih black hole in near 
and inner zone coordinates; e n is the perturbation parameter used in the approximate 
solution g$ in C„; is the mass of the ^4th black hole; b is the orbital separation 
on E; A is the gravitational wavelength. The quantity b is usually defined as the black 
hole separation in the near zone, which we shall see coincides with the black hole 
separation in the inner zone up to uncontrolled remainders, namely 6 = 6 + Q(e\, £2). 
Technically, there is a fourth region beyond the near zone, but we neglect this here 
since it does not affect the study of transition functions. For a detailed description of 
this subdivision and a pictorial representation of these zones refer to [HI [2] • 
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These regions overlap clearly in two buffer zones, which on the t = i = slice 
can be defined by the following inequalities: O13 = C\ n C3 (mj <C ri <C 6) and 
O23 = C% H C3 (m 2 <C r 2 <C 6) , where is the radial distance from the nth black hole 
to a field point. Asymptotic matching and the transition functions act in these buffer 
zones, which can only be defined asymptotically and, thus, any statement regarding 
them must be interpreted in that sense. In particular, this implies that any quantity 
that is valid in the buffer zone need not be valid near the boundary of the buffer zone 
(i.e., as ta — ► b or rA — > f/14.) Also note that these buffer zones are actually 4- volumes 
and can be foliated by an infinite number of timelike junction hypersurfaces. 

Different approximations are used in each zone to solve the Einstein equation. 
In either inner zone, black hole perturbation theory allows us to obtain a tidally 
perturbed metric. Let us concentrate on inner zone 1 near black hole 1, since the 
metric in the other inner zone can be obtained via a symmetry transformation. In C\, 
the perturbed metric is given in isotropic corotating coordinates x a = {t, x, y, z} by 

0$ « H t + H sl n 2 (x 2 + f) + 2H st x^ (x 2 +f- z 2 ) , 

'y(6lz-6lx)+(x 2 -z 2 )5 2 , 

9^ *6 al H sl -H s2 ^, (39) 

where f2 is the angular velocity of the perturbation, e a 5 E is the standard Levi-Civita 
symbol with convention e^s = 1 and where 5? is the Kronecker delta. In (f3"9")) we used 
the shorthand 

/l-Mi/2f\ 2 / Afi\ 4 m 2 , /i\ 

H < =-{TTM^r) +2 V^) -W**(r)> (40) 

where Mi is the mass of the background black hole, m 2 is the mass of the binary 
companion that is causing the perturbation, m = M\ + m 2 is the total mass, 
f = (x 2 +y 2 + z 2 ) 1 / 2 and F 2 (-) is the second Legendre polynomial. Equation ((39)) then 
satisfies the linearized Einstein equations in C\. These equations are identical to (18) 
and (19) of paper I and solve the Einstein equations up to uncontrolled remainders of 
C(fi/6) 3 . We refer the reader to that paper for an explanation of the derivation of 
this metric. 

In the near zone, a post-Minkowskian expansion is used to find an approximate 
solution. This solution in corotating ADMTT coordinates x a = {t,x,y,z} [ (7)-(12) 
in paper I] is given by 

(41) 



(!) ~ it o b , U st 

9o-a ~ " H ^e- am x + — 
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where we introduced a post-Newtonian conformal factor 

mi m 2 

* = H H -, 

2n 2r 2 ' 

and where the post-Newtonian lapse and shift are given by 



(42) 



2- * 



+ 



m i 
n 

m 2 

^2 



(y\ - vi ■ nin\) - Av\ 
(i> 2 - t/ 2 • n 2 n\) - Av l 2 



eik3^x 



(43) 



In these equations, the radial distance to the Ath black hole is given by ta = 
[x\ + y 2 + x 2 ) 1 / 2 , where x\ = x — m 2 6/m and x 2 — x + mib/m. These equations 
solve the Einstein equations in C3 up to uncontrolled remainders of ©(m^/r^) 2 ', Once 
more, we refer the reader to that paper for an explanation of the derivation of such a 
metric. In these equations, the non-zero components of the velocities va and the unit 
vectors Ha are given in the t = slice by 



m 2 , 

= CO 0, 

m 



x 



tk 



mi 

= —UJ 0, 

m 

m 2 



TA 



where w is the post-Newtonian angular velocity given by 







(-" 3 ) 


m 


H 


-I 








Vm / 


~b 



(44) 



(45) 



with errors of 0{m/bf/ 2 . Equation (g5]) appears also in paper I and in (60) of |31| . 
where /i = mim 2 /m is the reduced mass of the system. 

Asymptotic matching was performed in paper I and a coordinate and parameter 
transformation was found such that (J7]) is satisfied. For buffer zone 1, O13, such a 
transformation to Ofm^/r^, (fA/b) 2 } is given by 

( m 2 



Ma = mA: 



9. 



x 

l ~~b 



m 2 



where rj^ is the Minkowski metric and 8^ v is the 4-Euclidean metric. 



(46) 
Note that 



to zeroth order as previously mentioned and, when describing certain 



figures, we will sometimes use them interchangeably. One might worry that ()46|) 
is not bounded above, since x a — > 00 as t — > 00. However, recall that the buffer 
zone is technically a 4-volume delimited by the boundaries of the regions of validity 
of the approximations. In this sense, time cannot go to infinity, because then the 
post-Newtonian metric would break down as t approaches the time of coalescence. 
Furthermore, when constructing initial data, t — and i — 0, and it is clear that 
t — > 00 is not allowed. In paper I, asymptotic matching is carried out to slightly 
higher order, but since the purpose of this section is to study transition functions 
we use (|46[) instead. The metric in inner zone 2 and the matching coordinate and 
parameter transformations in buffer zone 2, O23, can be obtained by applying the 
following transformation to (|39jl and (j46| : 

1^2, x->-x, y^-y, z~>z. (47) 
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In order to study these transition functions, we pick a particular physical system. 
We choose a system of equal mass black holes mi = m 2 = 0.5m separated by an 
orbital distance b = 20m. The black holes are non-spinning and located on the £-axis 
at x = 10m and x = —10m for the first and second black hole respectively. The black 
holes are orbiting in the counter-clockwise direction about the z-axis. The buffer zones 
are then given by Oa3 ■ 20 3> r^/m ^> 0.5. This system is used to study different 
transition functions, but we should note that the asymptotically matched metric of 
paper I is valid for a wide range of systems. 

We here keep the physical system fixed and pick several different transition 
functions to investigate how these functions change the satisfaction of the Einstein 
equations. We consider the following type of transition functions: 



Transition function 1 (TF 1) 



r o 



, r < r_ 



h(r) = 



tanh 



s ( , 2w/k , . . 
- tan(—— (r-r_)) 

7T V bw 



tan(2^r(r-r_)) 



(48) 



< r < r + 



,r>r + , 

where r_ and r+ determine where the transition begins and ends, w determines 
the size of the region, and s and q determine the slope of the transition function 
roughly when fx sw 1/2. This function is similar to that used in paper I and we 
use similar parameters, i.e. 



r_ = m, r + = b — m, 
w = 9m, q — 0.2, s — b/m. 
• Transition function 2 (TF 2) and 3 (TF 3) 

/ b r ~ cT 



/ 2 (r) = exp 



1 



exp 



b r — d 
2m w 



(49) 



(50) 



\2m w 

where w is approximately the size of the transition window and d is approximately 
the distance from the Ath black hole at which the derivative of the transition 
function peaks. For TF 2 we choose the following parameters 

w = 7m, d = 6m, (51) 

while for TF 3 we choose 



d = 8m. 



(52) 



For all these transition functions, the transition window w is defined approximately 
as the size of the region inside which the function ranges from 0.01 to 0.99. 

The purpose of this paper is not to perform a systematic study of the properties of 
transition functions, but to illustrate the theorems discussed in earlier sections with a 
practical example. In fact, as discussed in the introduction, transition functions were 
first introduced in [51 [S] , but the properties and conditions that these functions must 
satisfy were not explored. In this paper, we are in essence following up on previous 
work and providing further details that might be relevant to future investigations. In 
particular, we shall investigate how the gluing of asymptotically matched approximate 
solutions breaks down if improper transition functions are chosen. 
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Figure 1. This figure shows the first (TF 1), second (TF 2) and third (TF 3) 
transition functions as a dashed, dotted and solid lines respectively. Observe that 
these functions transition at different radii and at different speeds. 



The transition functions that we have chosen clearly have different properties. 
First, notice that the region where functions arc significantly different from zero or 
unity varies, because the size of their respective transition windows w is different. 
In particular, TF 3 has the smallest transition window, followed by TF 2 and then 
TF 1. Second, notice that the transition functions become roughly equal to 1/2 at 
different radii: r ~ 7m for TF 1, r « 6m for TF 2 and r « 8m for TF 3. Third, 
notice that the speed at which the transition functions tend to zero and unity is also 
different. In particular, note that TF 3 tends to zero the fastest, followed by TF 1 and 
finally TF 3. This speed is important in the definition of proper transition functions in 
(|14| , because it prevents contamination of possibly divergent uncontrolled remainders 
outside the buffer zone. Finally, notice that not all functions are analytic, since TF 1 
does not have a Taylor expansion about r = r± (recall, however, that analyticity was 
not required.) 

These different features of the transition functions can be observed in Figures Q] 
and [2] Figure Q] shows the behaviour of the different transition functions as a function 
of radius, while figure [2] shows their first and second derivatives. Observe in figure [T] 
that the transition occurs at different radii and that they transition at different speeds. 
In figure [2] we can observe the difference in the transitioning speed better. We can 
clearly see in this figure that the derivatives of TF 1 are the smallest, followed by TF 2 
and then TF 3. However, both TF 1 and TF 2 have derivatives that are consistently 
of order much less than 0(1), while TF 3 has first derivatives of 0(1) or larger. The 
inset in this figure shows the large derivatives of TF 3. Also note that the second 
derivatives of the first and second transition functions are consistently smaller than 
their first derivatives. This is a consequence of the size of the transition window. If 
such a window were chosen to be smaller, as in the case of TF 3, then the second 
derivative would become larger than the first. 
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Figure 2. This figure shows the first and second derivatives of the transition 
functions. The first derivative of the first [(TF l),r], second [(TF 2), r ] and third 
[(TF 3) >r ] transition functions are denoted with a dashed, dotted and solid black 
line respectively. The second derivative of the first [(TF X) rT ], second [(TF 2) irr ] 
and third [(TF 3). rr ] transition functions are denoted with a dashed, dotted and 
solid cyan (or light gray on black and white print) line respectively. In the inset 
we zoomed out to show better the derivatives of the third transition function. 



With these transition functions we can now construct pure joined solutions via 
9^ = /(fi)/(f 2 )<$ + [1 - /(n)] g$ + [1 - /(f a )] gf}, (53) 

(12) 

where we assume g^ has been transformed with the coordinate and parameter 
transformation of (f46|) and where /(•) can be any of TF 1, TF 2 or TF 3. Equation 
(|53| is an extension of (fTTj) applicable to two buffer zones. Extensions to more than 2 
buffer zones are also straightforward. 

Whether the joined metric of ([53")) satisfies the Einstein equations to the same 
order as the approximate solution depends on the transition function used. We expect 
the metric constructed with TF 1 to generate small violations because it is a proper 
transition function that clearly satisfies the differentiability conditions required in 
Theorems 1 and 2. On the other hand, TF 2 and 3 do not satisfy these conditions 
because TF 2 is not a proper transition function and TF 3 violates condition (m) and 
(Hi) of Theorem 1. This behaviour, however, is not clear by simply looking at metric 
components. In figure[3]we plot the determinant of the spatial metric along the x-axis, 
corresponding to the x harmonic (near zone) coordinate, with t = 0, y = and 2 = 0. 
This axis is where joined metrics are glued together with different transition functions. 
This determinant gives a measure of the volume element on X and, as one can see 
in the figure, the difference when different transition function are used is small. This 
behaviour is not unique to the x-axis, but is actually observed along the other axis as 
well. When this is the case, we only show the behaviour along the x-axis in order to 
avoid redundancy. 

Although the volume element computed with different transition functions is 
similar, metrics joined with different transition functions do not satisfy the Einstein 
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Figure 3. This figure shows the determinant of the spatial metric along the x- 
axis with TF 1 (dashed line), TF 2 (dotted line) and TF 3 (solid line.) Observe 
that the differences in the global metrics are small if at all visible. In the inset 
we zoom to the region where the curves look the most different. 



equations to the same order. In figure HJ we plot the 4-Ricci scalar calculated with the 
global metrics joined with different transition functions. This plot is representative of 
the behaviour of the 4-Ricci scalar in the entire domain, although here we plot it only 
along the x-axis. We also plot only the region x/m > 10, to the right of where BH 1 
is located at t = 0, because the behaviour of the 4-Ricci scalar is symmetric about 
x — and the transition functions are symmetric about Ta = 0. Thus, the behaviour 
of the 4-Ricci scalar in other regions of the domain is similar to that shown in figure [4] 
(see paper I for contour plots of some of these quantities.) 

There are several features in figure |4] that we should comment on. First, 
observe that the is everywhere smaller than the uncontrolled remainder in the 
approximate solutions (roughly 0(m/b) 2 m 0.0025 in the buffer zone for the system 
considered) when the joined metric is constructed with TF 1. However, this is not the 
case for the metrics constructed with TF 2 and 3. For those metrics, ^'R has spikes 
close to ta — (x w 10m) and ta = 8m (x « 18m) for TF 2 and 3 respectively. 
The spike resulting from the metric constructed with TF 3 is associated with the 
small size of its transition window, which forces large derivatives in the transition 
function, thus violating conditions (ii) and (Hi) of Theorems 1 and 2. The spike in 
the metric built with TF 2 is related to the fact that this function does not tend to 
zero faster than the uncontrolled remainders of the post-Newtonian near zone metric 
close to BH A. As we discussed earlier, the rate at which transition functions tend 
to zero and unity is important to avoid contamination from divergent uncontrolled 
remainders in the approximations. In this case, TF 2 is of 0(m/b) 3 as r/ra — > m/b, 
while the post-Newtonian Ricci scalar diverges as m/r 3 as r — > and is in fact of 
0(b/m) 3 as r/m — > m/b. Since the transition function is not able to eliminate this 
contamination from the near zone metric, the Einstein equations are violated near 
either black hole for the metric constructed with TF 2. From a physical standpoint, 



Frankenstein's Glue: 
0.5 r 



0.4 



0.3 



(4) 



R0.2 



0.1 



Transition functions for approximate solutions 



22 



0.002 
0.001 


-0.001 
-0.002 



n 1 r 



/ s. 



i 1 r 



T 1 
T2 
T3 



_J I I I I L 

12 14 16 1 



1 I I I I I I I I I I I I I I I I 3 I I I I I I L 

10 11 12 13 14 15 16 17 18 19 20 21 22 

x/m 



Figure 4. 4-Ricci scalar along the x-axis for the global metrics joined with TF 1 
(dashed line), TF 2 (dotted line) and TF 3 (solid line.) Observe that the 4-Ricci 
scalar is small everywhere for the metric constructed with TF 1, but it has spikes 
for that built with TF 2 and 3. The inset zooms to the region x £ [11, 22] so that 
the difference between the 4-Ricci scalar calculated with TF 1 and 2 are more 
noticeable. 



the failure to use appropriate transition functions to create a pure joined solution 
introduces a modification in the matter-sector of the Einstein equations, namely the 
artificial creation of a shell of matter. One can straightforwardly see from (O and (|32| 
that the stress-energy tensor of this shell depends on the non-vanishing derivatives of 
the transition function. 

A deeper analysis of the spike in ^ R due to the joined metric with TF 2 reveals 
that the conditions found are sufficient but formally not necessary. The 4-Ricci scalar 
constructed with the post-Newtonian metric (^'Rpn) diverges close to BH A, since 
this approximation breaks down in that region. Therefore, if the transition function 
does not vanish identically, or faster than the divergence in 

<- 4) R PN , then ^R will 

present a spike. This spike, however, would disappear if the transition function decayed 
to zero faster than the divergence in ^Rp^. In this sense, eventhough the conditions 
discussed here are not necessary since the definition of a proper transition function 
could be weakened, they are certainly sufficient and universal. 

The second feature we should comment on is the behaviour of ^R in the region 
x/m £ [11,17]. As one can observe in the inset, the transition functions indeed 
introduce some error in this region. However, note that this error for the metric 
constructed with TF 1 is smaller than the uncontrolled remainders in the approximate 
solutions. Also note that this error is independent of the location of the center of the 
transition window. In other words, asymptotic matching is performed inside of a 
buffer zone and not on a patching surface, which means that approximate solutions 
can be glued with transition functions anywhere in the buffer zone away from the 
boundaries and not just at a specific 2-surface. As one can see from the figure, the 
metric constructed with TF 1 has the smallest error in this region, followed by that 
built with TF 2 and 3. This fact is not too surprising because, as shown in the proof 
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of Theorem 1, the error introduced by the transition functions in the calculation of 
the Ricci tensor scales with the first and second derivatives of these functions. Finally, 
note that the error introduced by these functions seems to be correlated to both the 
size and functional form of the second derivative of the transition functions, as one 
can see by comparing the inset of figure [2] to figure [4j The reason for such similarity, 
however, is beyond the scope of this paper. 

Once an approximate global 4-metric has been found, we can use it to construct 
initial data. These data could be given for example by 

h ab = f(h)f(f 2 )h$ + [1 - /(fx)] h% + [1 - /(fa)] 

K ab = f(ft)f(f 2 )K^ + [1 - /(n)] K$ + [1 - /(fa)] tf£>. (54) 

Equation ([54]) is simply a generalization of ([31]) for 2 buffer zones. The extrinsic 
curvatures for the inner zones and near zone are provided explicitly in paper I. 

Are we allowed to construct data with transition functions in this way? There 
might be some doubt as to whether this is valid, since the fact that the joined 
metric satisfies the Einstein equations does not necessarily guarantee that we can 
use transition function to construct the data itself. In particular, there might be some 
worry that the extrinsic curvature calculated from ([53]) generally contains derivatives 
of the transition functions that ([5l]) neglects. Theorem 2, however, ensures that 
this construction is indeed valid, provided the transition functions satisfy the same 
differentiability conditions proposed in Theorem 1. This is because the derivatives of 
these functions are then small and, in particular, the terms that are proportional to 
them are of the same order as the uncontrolled remainders in the approximations. 

These expectations can be verified in figure [5] where we plot the xy-component 
of the extrinsic curvature along the x-axis constructed both via ([54]) (referred to as 
glued T 1, 2 and 3 in the figure) and by direct differentiation of ([53]) (referred to 
as full T 1, 2 and 3 in the figure.) We plot only the cry-component because this is 
the dominant term of this tensor along the £-axis for the system considered and it 
shows the main differences in using different transition functions. We have checked 
that other quantities, like the trace of the extrinsic curvature, behaves similarly. For 
contour plots of the extrinsic curvature refer to paper I. 

As in the case of the Ricci scalar, there are several features of figure O that we 
should discuss. First, observe that that the extrinsic curvatures constructed with 
TF 1 agree in the buffer zone up to uncontrolled remainders. These remainders are 
of C(m/6) 3 / 2 « 0.01 in the buffer zone, because here we have used a coordinate 
transformation from the matching scheme that is valid only up to 0(m/b). The inset 
zooms to a region close to the outer boundary of the buffer zone in order to show this 
agreement better. The humps in this region are produced by the non-vanishing Lie 
derivatives of the transition functions. Observe that, as expected, these humps are 
smallest for TF 1, followed by TF 2 and 3. Our analysis suggests that if we performed 
matching to higher order, the agreement would be better and the size of the humps 
would decrease. The agreement is good for the curvatures constructed with TF 1, but 
the curvatures built with TF 2 and 3 have strong spikes roughly near the boundary 
of the buffer zone. These spikes arise because TF 2 violates the definition of a proper 
transition function, while TF 3 violates conditions (ii) and (Hi) of Theorem 2. 

We have then seen in this section how Theorems 1 and 2 can aid us in constructing 
transition functions. Even though we have not shown the error in the constraints, we 
have checked that this error presents the same features as those shown in figure 2] We 
have also seen the importance of restricting the family of transition functions to those 
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Figure 5. xj/-component of the extrinsic curvature along the x-axis, constructed 
via ((54J (black) for TF 1 (dashed line), TF 2 (dotted line) and TF 3 (solid line.) 
We also plot this component constructed by direct differentiation of H53H (cyan, 
or light gray on black and white print) for TF 1 (dashed line), TF 2 (dotted line) 
and TF 3 (solid line.) The inset zooms to the region x £ {15,22} so that the 
differences are more noticeable. 



that satisfy conditions (m) and (iii) of Theorems 1 and 2, as well as the definition of 
a proper transition function. The construction of pure, or mixed, joined solutions can 
then be carried out with ease as long as transitions functions are chosen that satisfy 
the conditions suggested here. 

5. Conclusion 

We studied the construction of joined solutions via transition functions. In particular, 
we focused on pure joined solutions, constructed by gluing analytical approximate 
solutions inside some buffer zone where they were both valid. The gluing process 
was accomplished via a weighted-linear combination of approximate solutions with 
certain transition functions. We constrained the family of allowed transition functions 
by imposing certain sufficient conditions that guarantee that the pure joined solution 
satisfies the Einstein equations to the same order as the approximations. With these 
conditions, we formulated and proved a theorem that ensures that the joined solution 
is indeed an approximate solution to the Einstein equations. We extended these 
conditions to projections of the joined solution onto a Cauchy hypersurface. We 
verified that the data on this hypersurface can itself be constructed as a weighted- 
linear combination with the same transition functions as those used for the joined 
solution. We proved that if these functions satisfy the same sufficient conditions, the 
data is guaranteed to solve the constraints of the Einstein equations to the same order 
as the approximations. 

We explicitly verified these theorems numerically by considering a binary system 
of non-spinning black holes. The approximate solutions used were a post-Newtonian 
expansion in the far field and a perturbed Schwarzschild solution close to the black 
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holes. We considered three different kinds of transition functions, two of which violated 
the conditions proposed in the theorems. The joined solutions constructed with these 
transition functions were shown to produce large violations to the Einstein equations. 
These violations were interpreted as the introduction of a matter-shell, whose stress- 
energy tensor was shown to be related to derivatives of the transition functions. The 
joined solution constructed with the transition function that did satisfy the conditions 
of the theorems was seen to introduce error comparable to that already contained in 
the approximate solutions. We further verified that projections of this joined solution 
to a Cauchy hypersurface can also be constructed as a weighted-linear combination of 
projections of approximate solutions with transition function. These joined projected 
solutions were seen to still satisfy the constraints of the theory to the same order as 
the approximate solutions provided the transition functions satisfied the conditions of 
the theorems. 

The glue proposed here to join approximate solutions has several applications to 
different areas of relativistic gravitation. Initial data for numerical simulations could, 
for example, be constructed once approximate solutions to the system have been found 
and asymptotically matched inside some buffer zone. In particular, initial data for a 
binary system of spinning black holes could be generated by asymptotically matching 
and gluing a tidally perturbed Kerr metric |32j to a post-Newtonian expansion (see 
[55].) One could also study the absorption of energy and angular momentum [3"4"j . 
as well as the motion of test particles [35], in the spacetime described by such pure 
joined solutions. Finally, one could extend the scheme developed here to formalize 
the construction of mixed joined solutions. Such solutions could, for example, be 
composed of a post-Newtonian metric glued to a numerically simulated metric or 
a semi-analytical approximate metric [101 1 1 1 L 1361 137j . Such mixed joined solutions 
could then be used to construct reliable waveform templates for extremely non-linear 
physical scenarios. 
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